Astronomy & Astrophysics manuscript no. cignoni2col 


© ESO 2008 


February 5, 2008 





The Galactic Halo density distribution from photometric survey 

data: results of a pilot study 

M. Cignoni 1 , V. Ripepi 1 , M. Marconi 1 , J. M. Alcala 1 , M. Capaccioli 1 , M. Pannella 2 , R. Silvotti 1 



O 

o 

(N 

> 
O 



> 
On 



O 

S3 

Of 
6 

CO 



1 Osservatorio Astronomico Di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy 

2 Max-Planck-Institut fiir extraterrestrische Physik, Giessenbachstrasse, Postfach 1312, D-85741 Garching, Germany 



X 



Received; accepted 



ABSTRACT 



Aims. Our goal is to recover the Galactic Halo spatial density by means of field stars. To this aim, we apply a new technique to the 
Capodimonte Deep Field (OACDF, Alcala et al. 2004), as a pilot study in view of the VLT Survey Telescope (VST) stellar projects. 
Considering the unique chance to collect deep and wide-field photometry with the VST, our method may represent a useful tool 
towards a definitive mapping of the Galactic Halo. 

Methods. In the framework of synthetic stellar populations, turn-off stars are used to reconstruct the spatial density. The determination 
of the space density is achieved by comparing the data with synthetic color-magnitude diagrams (CMDs). The only assumptions 
involve the IMF, age and metallicity of the synthetic halo population. Stars are randomly placed in the solid angle. The contributions 
of the various Monte Carlo distributions (with a step of 4 kpc) along the line of sight are simultaneously varied to reproduce the 
observed CMD. 

Results. Our result on the space density is consistent with a power-law exponent n x 3 over a range of Galactocentric distances from 
8 to 40 kpc. 

Key words. Galaxy: halo, Galaxy: structure, Galaxy: stellar content, Stars: Hertzsprung-Russell (HR) and C-M, Methods: statistical 



1. Introduction 

The astrophysical interest for the spatial distribution of the 
Galactic halo is fueled by many important questions concern- 
ing the origin of our Galaxy and its cosmological implications. 
In the classical scenario (Eggen, Lynden-Bell & Sandage 1962) 
the origin of the halo is identified in a rapid collapse of a proto- 
galactic cloud. An alternative scenario (Searle & Zinn 1978) sug- 
gests a different history: the Galaxy (and the halo) was built from 
smaller blocks, where the stars have already started developing. 
Signs of the accretion process are often identified in the observed 
streams of gas and stars, remnants of tidally-disrupted satellite 
galaxies. The emerging picture, suggested by recent results, has 
an hybrid nature: the outer halo may have been formed from the 
accretion of satellite galaxies, while the inner halo is thought to 
be a relic of a dissipative collapse. An evidence in favour of this 
scenario is that the axis ratio seems to vary from inner to outer 
regions (see e.g. Sommer-Larsen & Zhen 1990 Gould, Flynn & 
Bahcall 1998). Variations of the [a/Fe] ratio abundance are also 
observed: stars with lower [a/Fe] tend to stay at larger orbits 
than the richer ones, suggesting that the outer parts of the halo 
may have been accreted by intruder galaxies (see e.g. Nissen & 
Schuster [T997l 

In this context, the observation of a particular decline in 
the density slope could help to identify the right scenario. The 
mass density profile has been extensively studied, both using N- 
body simulations and matching suitable stellar tracers. Different 
power law indexes n (r~") are found: 

- n « 2 is required to explain the flat rotation curve; 



Send offprint requests to: M. Cignoni 



- n between 2.5 to 3.8 (e.g. Avila-Reese et al. 119991 ) is sup- 
ported by N-body simulations for outer halos; 

- n « 3 is found from satellite galaxies (SDSS data) orbiting 
around isolated galaxies (Prada et al. 120031 ): 

- n « 3.3 is found from local halo stars (Sommer-Larsen & 
Zhen ll99Tjb in the galactocentric interval 8-20 kpc; 

- n w 3 is derived from RR Lyrae stars (Vivas & Zinn 2003 
Ivezic et al. 2004]) in the galactocentric range 4-60 kpc; 

- n * 2.7 (Siegel et al.|20® and n ~ 2.5 (Robin et al. |2000b 
from star counts. 

These estimates suggest that different "candles" may yield 
different results. In particular it is hard to reconcile the isother- 
mal profile with the steeper declines derived both from simula- 
tions and from stellar tracers. The goal of this paper is to ex- 
tend these works, including also turn-off stars. In particular, we 
propose to recover the halo spatial density by comparing the ob- 
served CMD with a synthetic halo population (Monte Carlo gen- 
erated). Our result is a pure density distribution of stars along the 
line of sight (no analitycal a priori halo density is adopted). Only 
a posteriori fitting is performed on the recovered density. 

Section 2 is dedicated to the sample selection. In Section 3 
we discuss the methodology. In section 4 the reconstruction of 
the space density is fitted with a power law and the robustness of 
the result is tested against different assumptions on the adopted 
IMF and metallicity. Also the possibility of a metallicity gradient 
is considered. The conclusions close the paper. 



2. Data selection 

The OACDF is a multi-color imaging survey mainly for extra- 
galactic studies, acquired with the Wide Field Imager (WFI) at 



2 



M. Cignoni, V. Ripepi, M. Marconi, J. M. Alcala, M. Capaccioli , M. Pannella, R. Silvotti: Structure of the Galactic Halo 




FLUX_RADIUS 

Fig. 1. Visual apparent magnitude versus the sextractor param- 
eter FLUXJiADIUS (two different CLASS JSTAR selections are 
shown: filled circles are all objects, open circles are objects with 
stellarity larger than 0.9). 



the ESO 2.2-m telescope. This field is located at high galactic 
latitude (Z ~ 250, b ~ 50), making available a collection of 
halo stars with low contamination of disk and thick disk stars. 
The field of view covers an area of about 0.5 deg 2 in the B, 
V, R optical bands, reaching the following limiting magnitudes: 
Bab ~ 25.3, Vab ~ 24.8 and Rab ~ 25.1. The completeness 
magnitude limit is found at Vab ~ 24.0 (for a S/N ratio 10). 

2.1. Star/galaxy separation 

Contamination by galaxies dominates counts at faint magni- 
tudes. Thus, adequate star/galaxy separation is essential to in- 
fer the Galactic structure. In order to remove the bulk of 
extra-galactic contaminants, we applied a filtering based on 
the SExtractor (Bertin & Arnouts 1996) stellarity parame- 
ters (CLASSuSTAR) and half-light radius (FLUXJiADIUS). The 
point-like sources are identified by choosing objects with stel- 
larity larger than 0.9 (simultaneously in B and V filters). Beyond 
the limit V = 23, the photometric error in the V filter (see 
Alcala et al. 120041 1 overcomes 0.05 mag, making difficult a clean 
separation. This finding is evident in Figure Q] where the stel- 
larity selection is performed in a diagram V magnitude versus 
FLUXJiADIUS: selecting stars with high stellarity removes ob- 
jects with high radius (candidate galaxies). A side effect of this 
procedure is a loss of very faint stars (V > 23). Therefore, we 
limit our analysis to objects brighter than V = 23. 

The prescription on the stellarity (plus the cut in magnitude) 
does not ensure the complete removal of all no-stellar objects. 
Some point-like sources like quasars may be still present. To 
solve this problem we perform a cut in B-V (rejecting all objects 
bluer than 0.3) with a minimal effect on the stellar population 
(about 60 objects removed)^. 

2.2. Choice of a stellar tracer 

In order to obtain the spatial density of the Galactic halo one 
must identify stars in a precise evolutionary phase and use them 
to trace the observed CMD density. The strong requirement on 

1 Halo stars are old and metal poor: assuming the typical B-V color 
of old globular cluster, we don't expect to find stars bluer than B - V x 
0.3 - 0.4 (except a negligible fraction of white dwarfs and HB stars). 
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Fig. 2. Color-magnitude diagram of the point-like sources in the 
OACDF. The turn-off stars discussed in the text are represented 
with open symbols. 



the tracer used to reconstruct the space distribution is that it must 
be univocally linked to the distance, once metallicity, mass and 
age are fixed. In this context, red giants have an important neg- 
ative effect, because of the degeneracy (in color) with main se- 
quence stars. In contrast, turn-off stars have many advantages: 

- They can be univocally selected (definitely bluer than halo 
red giants); 

- They are readily distinguishable from the disk stars in the 
main sequence (only very young and massive disk stars may 
have similar colors, but they are rare); 

- They are numerous, in fact they outnumber HB stars and red 
giants by a factor > 100. 

On this basis, we decided to use turn-off stars as tracers of 
space distribution. In order to select turn-off stars we define all 
object in our sample with B-V color in the range 0.3-0.6 (see 
figure [2]). The final data contains about 600 stars. 

2.3. Explored distances 

In what follows we present results for stars in the range 18 < 
V < 23. This achieves two goals: the bright limit was set to 
avoid thick disk stars, while the faint limit minimizes the con- 
tamination of no-stellar sources. Moreover, correlations due to 
photometric errors are also reduced. These boundaries guarantee 
a theoretical range between 4 kpc and 60 kpcQ (heliocentric). 

3. The model 

A new method for calculating and displaying patterns of Galactic 
halo stellar density is presented. The recovered distribution 
N(r Q ) is the number of stars for heliocentric interval 5r & and 
solid angle u>. To recover the spatial density behind the observa- 
tional CMD we need to model a synthetic halo population with 
various theoretical ingredients. Then the match between theo- 
retical and observational CMDs is quantified (see Cignoni et al. 
2006 for details). 



2 A simple population with metallicity Z=0.001 and 12 Gyr, the tum- 
off luminosity is near to M v ~ 4 (lower metallicities lead to a brighter 
limit). In this case, a minimum luminosity V = 23 implies a maximum 
distance around 60 kpc. 
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Following a Monte Carlo technique, the theoretical CMD is 
populated with a large number of stars randomly built by means 
of a set of theoretical ingredients. Ignoring chemical evolution, 
the synthetic halo population is built near coeval and mono 
metallic. The adopted star formation rate is constant between 10 
and 12 Gyr (age estimations for globular clusters indicate values 
between 1 1 Gyr, Chaboyer et al. 1998 and 13 Gyr, Hansen et al. 

2002). The metallicity is fixed at [Fe/H] 1.6 (see e.g. Ryan 

& Norris 1991) with a Gaussian spread (o-[f e /H] = 1.0). The 
adopted IMF is a power law with a Salpeter index. Once masses, 
ages and metallicity are known, colors and magnitudes are ob- 
tained by interpolating a library of stellar evolutionary tracks 
(Pisa Evolutionary Library, Cariulo et al. 120041 ). 

In order to reconstruct the halo structure, this synthetic pop- 
ulation is placed at different distances, randomly placing stars 
in the j-th heliocentric interval [c/ .,, d Q j + 4 kpc]. The syn- 
thetic diagrams are corrected for reddening and extinction ac- 
cording to Schlegel et al. 119981 maps (for the OACDF field, a 
E(B -V)~ 0.05 is adopted). The final product is a base of par- 
tial color magnitude diagrams, each one representing the same 
population but at different distance moduli. Transforming these 
CMDs in two dimensional histograms (bin size: 5V = 0.05, 
mag, 5{B - V) = 0.05 mag), we obtain a base in a 2-D space 
(V versus B - V) and we can write every CMD as a linear com- 
bination of these histograms. In numbers, a generic CMD is so 
built: 



ii = Yj a i Ci i 



(1) 



where ;«,■ is the number of star in the final CMD in bin i, aj 
is the stellar rate in the distance slab j, and cy is the number of 
stars in the bin i owning to the partial CMD /'. 

The combination of coefficients aj that minimize a Poisson 
based likelihood is searched through a simplex algorithm. The 
uncertainties on the final coefficients are found by using a boot- 
strap technique. The final result is the distribution of stars along 
the line of sight in the direction of the OACDF. However, real 
data are observed in solid angle, therefore our result is not a den- 
sity yet. After correcting for the solid angle, the space density as 
a function of the Galactocentric distance is found by converting 
from heliocentric distances to Galactocentric distances. 
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Fig. 3. The recovered density distribution (the density scale is 
somewhat arbitrary). Error bars are 1 cr. For comparison, differ- 
ent power law densities are also shown (with the labeled expo- 
nents). 
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Fig. 4. Reduced^- 2 versus the n exponent. 

interval for «, these results confirm a stellar halo with a decline 
steeper than an isothermal profile (which is the accepted profile 
for the dark halo). Indeed, values n lower than 2 give axl more 
than 4 times higher than the n = 3 result. 



4. Results 

The recovered spatial distribution is shown in Figure[3](symbols 
with error bars). We have to stress the exact nature of our re- 
sult. In fact, the OACDF field is not along a galactic radius, but 
it is rather oblique. Thus, comparing our result with simple ra- 
dial profiles means that we are assuming an axially symmetric 
Galaxy. Moreover, due to the large uncertainties of the further 
bins, our analysis is truncated at the galactocentric distance of 
40 kpc. 

A parametrized model of the mass distribution is fitted to the 
recovered density. The halo is described by a power law r~" (r 
galactocentric radius), with axis ratio equals to 1 (round halo). 

Different hypothesis on the power law index n (see figure[3j 
are explored through a chi square fitting technique. 

Figure [4] illustrates the reduced xl 0^ 2 /degrees of freedom) 
versus the power law index n: the minimum is placed near n = 
2.9 with;^ 2 « 1.0. Although it is difficult to define a confidence 



4.1. Alternative scenarios 

In order to avoid biases due to the assumed halo properties, 
namely the IMF, the metallicity Q and the axis ratio, we have 
repeated our recovery by using different scenarios. 

4.2. IMFs 

Figure [5] shows the effect of different IMFs. The result is quite 
independent of the precise IMF (dN/dm <x m~ s ) exponent s, 
probably due to the short mass spectrum involved in the halo 
population. The cuts at V — 23 and B — V = 0.6, combined 

3 The uncertainties on the recovered densities may hide correlations 
among bins (incompleteness, external contaminations and wrong metal- 
licitiy may produce systematic effects; see e.g. Cignoni et al. 2006), 
making the^ 2 probability quite useless. 

4 In principle also the assumed age is critical. However, the small 
masses involved (lower than 0.8 M ) imply a weak sensitivity to the 
age. 
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Fig. 5. Reduced^ 2 versus index n for different assumed IMF (the Fig. 7. Reduced y 2 versus index n for several hypothesis on the 
explored exponents are labeled). axis ratio (AR= 0.6, 0.8, 1.0). 
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Fig. 6. Reduced x 2 versus index « for several assumed halo 
metallicities. The metallicity distribution is always a Gaussian 
with the labeled mean and dispersion; the label "grad" indicates 
a model with a spatial gradient - see text. 



with the typical low metallicities involved in the halo, imply a 
detectable minimum mass of about 0.6 M . Then, the age of the 
halo, that is higher than 10 Gyr, implies a maximum mass (for 
stars still in main sequence) around 0.8 M . In conclusion, the 
total mass spectrum is only =s 0.2 solar masses: too short to dis- 
tinguish between different IMF slopes. Photometric errors, the 
presence of different metallicities and low number statistics can 
largely sweep out an IMF effect. 

4.3. Metallicity effect 

A variation in the metallicity spread has not a significant impact 
(figure|6]shows the result for A[Fe/H] = 0.3,0.5, 1.0). In fact, a 
metallicity spread introduces a blurring on the partial CMDs and 
this causes an additional random uncertainty on the recovered 
density. However, this effect is not systematic. 

The result is different when the mean metallicity of the 
model is changed. In particular, increasing the metallicity to 
[Fe/H] = -1.0 shifts the most probable n exponent to lower 
values. The new best value is close to n = 2, but the relative;^ 2 , 
is substantially higher than the previous estimates (^ 4.5). 

Finally, Figure [6] shows also the resulting xl when an ex- 
treme spatial gradient (0.2 dex each 4 kpc, [Fe/H] = in the 



Solar position) is adopted (labeled with "grad"). As expected, 
the adoption of a metallicity gradient has an intermediate effect: 
the minimum xl is reached for n as 2.3, yielding a xl value 
around 1.3, which is comparable with those resulting from the 
first models. 

4.4. Axis ratio 

In this section the hypothesis of a round halo is abondoned. 
Although the axis ratio (AR) is not significantly constrained by 
a single field, two alternative values are explored (AR = 0.8 and 
0.6). Figure |7] shows the results. Decreasing the AR determines 
a lower best value for n. The extreme choice AR=0.6 implies a 
value n « 2.2. It is noteworthy that the best n, independently of 
the precise AR, is caracterized by similar chi square values: for 
a single field, the degeneracy between axis ratio and power law 
index does not allow to disentangle the two effects. 



5. Conclusions 

We have combined a star count model with an updated stellar li- 
brary to reconstruct the Galactic Halo stellar density in the direc- 
tion of the Capodimonte Deep Field. Our results confirm a steep 
decline. In particular, when a metallicity [Fe/H] = -1.6 ± 1.0 
is adopted, a power-law index n « 3 gives a reasonable fit out 
to » 40 kpc from the galactic centre. This is quite different from 
the value n « 2, implied by the flat rotation curve. However, gas 
and turn off stars may trace different halo regions. First, the mea- 
surements of rotation velocity from gas are reliable up to about 
20 kpc from the Galactic Center while turn off stars reach 40 
kpc, suggesting that our exploration is deeper in the outer halo. 
Second, the OACDF field is high on the Galactic plane while the 
gas is locked on the plane. 

The value n » 3 is also larger than the Robin et al. 2000 and 
Siegel et al. [2002 findings, that analyzed several galactic fields 
(generally shallower than OACDF) using similar stellar tracers. 
Basically, our study confirms the results obtained with RR Lyrae 
stars. 

An higher metallicity ([Fe/H] = -1.0) can lower our best 
value to the isothermal n » 2, but the corresponding xl 1S sig- 
nificantly worse. Moreover, the presence of an extreme chemical 
gradient (xl values similar to the [Fe/H] — -1.6 + 1.0 scenario) 
cannot be ruled out. 
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A lower index n is also obtainable by adopting smaller axis 
ratios. In these cases, the minimum xl are close to case with 
round halo (AR=1.0), leaving the possibility of a lower index n 
(n = 2.2 for AR=0.6). 

This work represents a pilot program for a larger project de- 
voted to the explorat ion of the VST halo fields (STREGA® VST, 
see Marconi et al. 2006). The possibility to perform multi- 
direction comparisons will allow us to reconstruct the 3-D pat- 
tern of the halo and to disentangle the degeneracy between the 
axis ratio and the density power law slope. Moreover, the avail- 
ability of several photometric filters will provide a chance to im- 
prove both the comparison between synthetic and observational 
CMDs and the separation star/galaxies. 
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